A vortex dynamics approach to the Nernst effect in fluctuating superconductors 
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We present a new method to study the Nernst effect and diamagetism of an extreme type-II 
superconductor dominated by phase fluctuations. We work directly with vortex variables and our 
method allows us to tune vortex parameters (e.g., core energy and number of vortex species). We 
find that diamagnetic response and transverse thermoelectric conductivity (a^cy) persist well above 
the Kosterlitz-Thouless transition temperature, and become more pronounced as the vortex core 
energy is increased. However, they weaken as the number of internal vortex states are increased. 
We find that a^y closely tracks the magnetization {—M/T) over a wide range of parameters. 
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A number of experimental observations in superfluidity 
and superconductivity are best explained in terms of the 
statistical mechanics and dynamics of vortices. Exam- 
ples include the Kosterlitz-Thouless (KT) transition in 
superfluid films, and the "flux-flow" contribution to elec- 
trical resistivity in type-II superconductors. A natural 
extension of this paradigm involves the use of thermo- 
electric and thermal transport experiments as a probe 
of vortex dynamics. This is especially pressing given the 
Nernst effect experiments in the cuprates [ij . The Nernst 
effect is defined as the appearance of a steady-state elec- 
tric field (Ey) when a system is placed in a perpendicular 
thermal gradient (V^jT) and magnetic field (H^) under 
open circuit conditions. In the cuprate superconductors, 
a large Nernst signal and diamagnetic response persists 
well above Tc, the superconducting transition temper- 
ature, and are especially enhanced in the underdoped 
regime of hole-doped materials. 

To explain these results, it has been argued that vor- 
tices are responsible for the large Nernst signal [T]: as a 
vortex drifts down the temperature gradient, it generates 
a transverse electric field. A vortex description above Tc 
is useful when the superconducting transition represents 
the loss of macroscopic phase coherence, while the ampli- 
tude of the order parameter remains large below a higher 
'mean-field' temperature scale T^ ^^ ^ 7c Hi- This idea 
is most likely to hold true in the underdoped regime of the 
hole-doped cuprates, where the reduced superfluid stiff- 
ness enables phase fluctuations to suppress Tc well be- 
low the mean-field transition temperature. This "vortex- 
liquid regime" regime Tc < T <^ Tc'^ [3| over a range 
of fields can be intuitively understood in terms of a di- 
lute fluid of vortices, but so far this appealing picture 
0, m has not lead to a quantitative understanding of the 
Nernst effect. 

There are special challenges involved in constructing 
a vortex based theory of thermoelectric transport that is 
consistent with the basic principles of statistical mechan- 
ics. In contrast to thermodynamics, and even electrical 



transport, where a theory of vortices interacting via a 
long range potential can be used, a thermal or Nernst 
transport calculation requires a purely local formulation 
where no such long range forces are explicitly present. In 
this Letter, we present such a local formulation, and use 
it to study thermoelectric transport directly in the vor- 
tex language. An advantage of our method is that vor- 
tex parameters such as core energy, or number of vortex 
species, can be tuned independently of other properties, 
and their impact on the Nernst signal and magnetization 
can be systematically studied. Moreover, it is possible to 
use our method to study the Coulomb gas efficiently even 
in the high density limit, where the number of pairwise 
Coulomb interactions is large. 

Previous theoretical work on the Nernst effect in 
cuprates 0, 0, H, S has studied the time-dependent 
Ginzburg-Landau theory, which includes fluctuations in 
the amplitude of the order parameter. In earlier work, 
three of us have studied the Nernst effect in an XY 
model fl3| , which should be relevant for the underdoped 
cuprates. Our present results approaching from the vor- 
tex viewpoint are consistent with this earlier approach. 
Remarkably, we find that the close quantitative connec- 
tion between the Nernst effect and diamagnetism ob- 
served in that and other papers, holds in the present 
study too, even when vortex properties are drastically 
modified. The effect of vortex core energy on diamag- 
netism and the transition temperature was studied in 
refs. 
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Method: Consider a two dimensional superconductor 
in the extreme type-II limit, in which the supercurrents 
are too feeble to modify the externally imposed field. 
Since we are interested in the 'vortex liquid' regime, we 
restrict the order parameter to live on the sites n of a 
square lattice and assume that all fluctuations arise from 
the phase: = |V'o|e'^"^*^ We map to the vortex rep- 

resentation by defining a dual electric field e = (ef,e^) 
on the bonds {i, x); + y) of a dual lattice, orthog- 
onal to the original lattice bonds, via: e = V9 x z where 
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V is a lattice derivative. The local field e is related to 
the local supercurrents via J = x e, where is the 
bare superfluid density. Vortices may live at the sites i of 
this dual lattice, and the integer vortex charge satisfies 
Gauss' law V • e = 27rni. We thus obtain a completely lo- 
cal Hamiltonian for the vortex fluid that retains all phase 
and vortex degrees of freedom: 
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(1) 



where Sc is the vortex core-energy and we use units where 
= l/27r. When the transverse electric fields are inte- 
grated out, the model reduces to the static 2D Coulomb 
gas. However, we stress that one must include both lo- 
cal supercurrents and vortices when dealing with thermal 
transport, so as to define a local energy density. The in- 
teraction between the vortices is mediated by the super- 
currents, and this maps simply to the interaction between 
the charges being mediated by the dual electric field in 
the Coulomb gas. 

The model is given a Monte Carlo dynamics that cap- 
tures the effect of random thermal fluctuations. Two 
distinct types of updates, charge and curl updates, cor- 
responding to the two degrees of freedom (vortices and 
supercurrents ) are introduced (Fig. [1]). During a curl 
update, a plaquette is chosen at random and a random 
electric curl is added to it. Such an electric flux con- 
figuration is purely transverse and is not accompanied 
by vortex creation. During a charge update, a lattice 
bond is chosen at random, and a vortex/ anti- vortex pair 
is added on the two sites connected to this bond. The 
electric fluxes are updated locally near the charges to sat- 
isfy Gauss' law. Note, this charge update may result in 
the motion a pre-existing vortex. In the simulation, each 
move is accepted with probabihty 1/ [1 + exp{AU/T)], 
where AU is the change in energy associated with the 
move and T is the local temperature at the center of the 
plaquette or bond for that move. By varying the rela- 
tive frequency of each type of trial, we have control over 
Dpfi, the phase diffusivity, relative to the vortex diffu- 
sivity Dy. In what follows, we work in the physically 
reasonable limit Dph S> Dy. After an attempt is made 
to update each bond and plaquette with such moves, a 
unit of Monte Carlo time passes in the simulation. A 
related method (without pair creation) has been used in 
Ref. to study charged polymers. 

When a charge move is attempted, the simplest way to 
satisfy Gauss' law is to add an electric flux e = 2TTf to 
the link r connecting the positive to the negative charge. 
However, such single bond updates usually produce large 
AU and thus are rarely accepted, resulting in too little 
vortex motion at low temperatures; hence an update that 
spreads the electric flux over several bonds is used. A 
simple example of such an update is shown in Fig.[T] The 
added flux is made curl- free and the move we actually use 



FIG. 1; Charge updates (left) and curl updates (right) sample 
the longitudinal and transverse degrees of freedom, respec- 
tively, of the dual electric field. 




FIG. 2: Dielectric response of the dual Coulomb gas, com- 
puted using Eq. [2] and (inset) rescaled for finite-size effects 
We use an L X L torus with L = 8, 10, 12, 14, 40, and 
core energy ec — 0.125. The transition temperature is the 
location where the rescaled data intersects the line y — AT 
(dotted line). 



involves a patch of 12 plaquettes that is one plaquette 
larger in all directions than indicated in Fig. [TJ 

Thermodynamics: We have studied the KT transition 
of the neutral 2D Coulomb gas by tracking the dielectric 
response function 
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At T = TxT, the dielectric response function satisfies 
the universal jump criterion [18]: 



£-1 = 
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T — Tj^rp 

T — T+ 



(3) 



Fig. [2] shows the results of the helicity modulus cal- 
culation using various system sizes and a core energy 
Ec — 0.125. We use finite-size scaling [l3| to identify 
Tkt = 0.22; the naive free-energy estimate is T^j, = 0.25 
in our units. In the remainder of this letter, we cite all 
temperatures in units of T^rp. 

Diamagnetism: When a magnetic field is applied, the 
resulting imbalance of vorticity leads us to consider a 
plasma of vortices in a static, neutralizing background 
(charge density rib = —B/^q), where $o = Stt is the flux 
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quantum in our units. To study diamagnetism we must 
permit vortices to enter and leave the sample, so we em- 
ploy a cylindrical geometry with open boundaries. Cur- 
rent flow near the boundaries is measured, which arises 
due to a surface depletion of vorticity. This has a lower 
free energy than the perfectly neutral system, and leads 
to diamagnetism. Deep inside the cylinder, beyond some 
distance xq from the edge, the supercurrents vanish in 
equilibrium. The magnetization can be obtained by in- 
verting the relation J = V x M on a cylinder whose 
axis is along the x direction. For physical clarity we 
present continuum formulas below, these can be read- 
ily transcribed into the appropriate lattice versions. We 
have: 



M 



dx (jy{x,y)) ^ ps 



dx (cx) 



(4) 



Thus, the magnetization is directly proportional to the 
work function of the dual Coulomb gas (energy cost of 
removing a vortex from the bulk of the system). Using 
Gauss' law we can also obtain: 



M = 27rps 



dx X i{n{x)) -B/^o) 



(5) 



Thus, the magnetization is also the total edge polariza- 
tion (dipole moment) per unit length. The polarization 
fields are non-zero only in the charge depletion region 
near the cylinder's edge. 

Nernst Effect: To determine the Nernst effect, we again 
make use of cylindrical geometry, and apply a tempera- 
ture gradient along the cylinder axis. In our simulations, 
we compute the transverse thermoelectric conductivity 



defined via {Jy 



-ctxy{—^xT). The thermoelec- 



tric conductivity is closely related to the Nernst signal: 
cnxyl <^xxi where Oxx is the electrical conductivity and we 
have made the approximation of vanishing Hall angle. In 
our calculations, axy is obtained by measuring the dual 
electric field e^, and using the relation axy = {e^) j'^ xT . 
Thus while oixy is an off-diagonal transport coefficient 
when written in terms of electrical currents, it is a di- 
agonal response function in the vortex representation: it 
is simply the vortex thermopower. We have verified that 
the net vortex motion vanishes once the steady-state in 
the thermal gradient is reached, as is required in a ther- 
mopower measurement. The quantity axy has the advan- 
tage that, unlike the Nernst signal bat, it does not have 
any explicit dependence on tmc, the Monte-Carlo time 
step, as can be seen from the Kubo formula for axy and 
from dimensional analysis 11 1 . 

Figure [3] shows the simulation results for axy and M 
(inset b) for apphed fields up to Bq = $o/(27ra^), where a 
is the lattice spacing, comparable to the zero temperature 
coherence length ^p- Both axy and —M/T are expressed 
in units of the 2D "quantum of thermoelectric conduc- 
tance," 2efcs /h. We show data for T > T^j,; below T^j,, 




FIG. 3: axy for a system of single-flavor (solid lines) and 
2-fiavor (dotted lines) vortices on a 20x15 cylinder, with core 
energy Ec ~ 0.5r]J-y. Temperatures shown are 1.0 (topmost 
A curves), 1.25, 1.5, and 2.0 T^^j. (bottom 2 curves). Inset 
a: Core energy dependence of a^y (solid lines) and —M/T 
(dashed lines) in units of 2e/cs//i for single-flavor vortices, at 
T = 2.5T^T, B = 0.7Bo showing a marked increase of both 
quantities with core energy. Inset b: Diamagnetism —M/T 
curves for the same temperatures as the main figure. (Units 
of 2efefl//i are used on the vertical axes of both insets.) 



our simulations encounter difficulties due to impaired 
vortex mobility, and furthermore, it is difficult to apply 
thermal gradients small enough to remain within linear 
response. Fortunately, T > T^rp is a regime of interest, 
since the Nernst signal persists well above Tc in the ex- 
periments of Wang, et al. [1]. The results presented here 
are in quantitative agreement with earlier computations 
involving the 2D XY model with Langevin dynamics ll| : 
in particular, both models have the feature that in the 
small magnetic field limit, a) axy and M diverge loga- 
rithmically as i? ^ for T < Tkt and 6) they increase 
linearly with B at small B when T > Tkt UM- Simi- 
lar features have been seen in experiments of Wang, et al. 

Moreover, both axy and M are detectable at temper- 
atures as high as 2Tkt in this vortex model, which is also 
consistent with 2D XY model results [ll]. At very high 
temperatures T » Tkt, the 2D XY model predicts that 
axy and M/T decay sharply as a power law in temper- 
ature. In the vortex model, however, the magnetization 
decreases even more rapidly at such high temperatures: a 
calculation based on the dual solid-on-solid model, shows 
that the magnetization of single-flavor vortices decays ex- 
ponentially as: M -2Ts\n{B/Bo)e-'^'^/p' We have 
checked that our numerical results agree with this expres- 
sion. Although we have not succeeded in finding similar 
expressions for axy, our numerical results indicate that 
axy also decays in this fashion at high temperatures and 
closely tracks the diamagnetism. 

Vortex Core Energy Dependence: The core energy de- 
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pendence of axy and diamagnetism —M/T are shown in 
Fig. [3] (inset a) at T = 2.5T^j,. So long as > T, 
both are found to increase with ej,. At this temperature, 
ttxy and —M/T track each other closely. With increas- 
and —M/T rise from near zero at 



— IQT'^j,, showing that the core 



mg core energy, a^y 
Ec = to 0(1) at e, 
energy has a dramatic impact on both of these quanti- 
ties in this regime. The dominant effect of the vortex 
core energy ec is that it enhances local superconducting 
correlations at short distances by increasing the cost of 
vortex fluctuations. The core energy enters directly in 
setting the "work function" for removing a vortex from 
the system, which is proportional to the magnetization. 
For tc ^ T, we thus expect —M oc tc, since more vor- 
tices are expelled near the boundaries in this limit. We 
have observed this in our simulations. The vortex- free 
boundary layer grows with tc and when it becomes com- 
parable to the thickness of the sample this causes strong 
finite-size effects. Thus we have not been able to reliably 
determine the bulk behavior for large values of the vor- 
tex core energy tc- We expect axy to saturate at a finite 
value in this limit - when all thermally-generated vortex 
fluctuations are suppressed, the remaining field-induced 
vortices respond to the thermal gradient in a way that is 
independent of the magnitude of tc- Although we do see 
some indications of the saturation of axy at large values 
of Ec, we have not been able to access this regime reliably 
due to the large finite-size effects on axy when tc » T. 

Dependence on the Number of Vortex Flavors: Several 
theories of cuprates predict additional degrees of free- 
dom associated with vortices, that endows them with 
a flavor index [i^j 21 1. Such an extension is readily 
incorporated in our formalism, the Gauss law is now 
V • e = 27r X]^=i where a is the flavor index . Our 
data shows a systematic dependence on the number of 
internal flavors of vortices. We find that both \M\ and 
axy of single-flavor vortices (solid curves in Fig. [3|) are 
systematically larger than that of vortices with 2 internal 
flavors (dotted lines). With multiple vortex species, there 
is additional conflgurational entropy associated with the 
internal degree of freedom. Therefore, the free energy 
cost of introducing a vortex into the sample is lower, and 
the system becomes less diamagnetic. It is more diffi- 
cult to understand, however, why axy decreases as the 
number of species is increased. In phenomenological dis- 
cussions, the transport coefficient axy is usually identi- 



fied with the vortex entropy 22|. Therefore, one may 



naively expect axy to increase with the number of in- 
ternal vortex fiavors. Moreover, axy is equivalent to the 
vortex thermopower, and if we neglected the logarithmic 
interactions between vortices, a simple calculation shows 
that the thermopower of an ideal gas of vortices is sim- 
ply the entropy per particle [23|] ; which increases linearly 
with the number of internal states associated with such 
particles. Our numerical results however are in sharp 
contrast to such expectations that axy ought to be pro- 



portional to an effective entropy per vortex. Instead, our 
results which have the opposite trend point towards a 
scenario that is far more complicated than an ideal gas 
approximation; interactions between vortices and ther- 
mally generated vortex-antivortex pairs invalidate such 
simple relations between axy and vortex entropy. 

Perhaps the most striking feature of our data is the 
fact that axy and —M/T closely track each other. In 
particular, at high temperatures, our results obey the 



relation 



-cM/T, where c « 1, and is difficult 



to determine accurately due to noise at large tempera- 
tures. Recently, it was shown analytically in Ref. [llj 
that for the 2D XY Hamiltonian with overdamped dy- 
namics, c = 1/2 at high temperatures. Furthermore, one 
obtains the same relation for Gaussian superconducting 
fluctuations [l^]- Here, we have shown that this relation 
is even more robust, since it hardly depends on the varia- 
tion of the core energy and the number of distinct vortex 
species so long as Cc is not much greater than T. This 
relationship which connects a transport coefficient and a 
thermodynamic quantity is indicative of a deep underly- 
ing (and as yet unknown) principle, that seems to hold 
true over a broad range of parameters. 

In conclusion, we have presented a new local method 
to study the thermodynamic and transport properties of 
2D vortex liquids. This method enables us to directly ob- 
serve the effect of tuning vortex parameters on axy and 
M . As an application of our method, we have systemat- 
ically studied the effect of modifying the vortex core en- 
ergy and the number of vortex species on diamagnetism 
and Nernst effect. In all cases, we have found that both 
quantities persist well above Tkt and axy closely tracks 
-M/T when Tkt < T ^ e^. We have provided the 
first detailed analysis of the Nernst effect of a vortex liq- 
uid that deals directly with vortex variables. The method 
presented here can be generalized to study thermal trans- 
port of vortex liquids in 3D, and can be extended to study 
quantum phase fiuctuations at low temperatures. 
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